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ABSTRACT 



The extreme speed and altitude ranges of which modern, high perfor- 
mance, supersonic airplanes are capable have become manifested in wide 
variations in airframe transfer function parameters as flight conditions 
change* Airframe control systems and stability augmentation devices 
utilized in these aircraft require some form of gain changing to pro- 
vide satisfactory stability and control characteristic over the complete 
flight envelope. In present day operational aircraft, system gains are 
scheduled as functions of airspeed and altitude, however these mechaniza- 
tions have seme shortcomings which are pointed out in this paper, and also, 
the desirability of self-adaptive aircraft control systems is discussed. 

A survey of adaptive techniques which have been applied to the airplane 
control problem is presented, and a self-adaptive pitch-stability aug- 
mentor which utilizes a digital computer in the adaptation loop is pro- 
posed, The proposed self-adaptive system is tested in conjunction with 
a digital three degree of freedom flight simulation of the pitch dynamics 
of a representative supersonic fighter-attack aircraft, and is found to be 
capable of performing the adaptation process in two cycles of airplane 
short-period oscillation. Requirements for practical mechanization in 
a flight vehicle and recommendations for future research are discussed. 
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- DTTROWCTIOil 

As the speed and c titude capabilities of modern high performance- 
supersonic aircraft have increased in recent years, the range of varia- 
tion of the dynamic characteristics of these aircraft have increased pro- 
portionately. A ten to one variation in damping ratio and natural fre- 
quency of some of the dominant modes of aircraft motion is not uncommon 
for modern airplanes with Mach 2 and 70,000 ft. altitude capabilities. 

As an example, for the airplane considered in this study, the damping 
ratio of the dominant longitudinal mode varies from £ = 0.445 at M = 0.9 
@ 10,000 ft. to£ = 0.064 at M = 1.2 <? 45,000 ft. 

Various stability augmentation devices have been used to decrease the 
range of variation of airplane parameters and to hold them as near as 
possible to optimum values. These devices, in their historical order of 
application, include fixed gain dampers, pilot-adjusted-gain dampers, 
speed-and-altitude- programmed -gain dampers, and se If -adapt ive- gain -adjust- 
ing dampers. A number of methods have been investigated for mechanizing 
self-adaptive stability augmentation systems for aircraft, utilizing vari- 
ous techniques which will be discussed later. And, the purpose of this 
study is to investigate the use of an airborne stored program digital 
computer for generation of the adaptation commands in a self-adaptive 
pitch damping system to be used in a high performance fighter- interceptor 
type aircraft. The dynamics of the airplane, including control system 
components, and the adaptation computer were simulated on a Control Data 
Corporation 1604 general purpose digital computer. 

Before beginning the discussion of adaptive control systems as applied 
to aircraft, it is felt that a bettor understanding of the system require- 
ments may be gained by a short review of the characteristic dynamic modes 
of a supersonic aircraft and the theory underlying the analysis and syn- 
thesis of these motions. 
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The airplane ikj i-j axis system which is generally used in airplane 
stability .raid c outre 3 studies is shown in figure 1, as three mutually 
perpendicular axes intersecting at the airplane center of gravity. The 
longitudinal, or X, axis lies in the airplane's vertical plane of sym- 
metry and is arbitrarily oriented by the airplane designer so as to de- 
fine a horizontal plane of symmetry of the fuselage. In almost all air- 
craft the X body axis does not coincide with the longitudinal principal 
axis (longitudinal axis about which the moment of inertia is least) due 
to the high location of the mass of the tail surfaces and due to the 
assymetric top-to-bottom configuration of the fuselage. Consequently, the 
longitudinal principal axis normally slopes downward from tail to nose and 
coincides with the X body axis only at the center of gravity. The dynamic 
unbalance of mass about the X axis thus causes the airplane to have an 
I product of inertia which will be important when cross coupling of air- 
plane motions is discussed. 

The Z body axis, or normal axis, is perpendicular to the X axis and to 
the wing plane and has positive direction downward. And, the Y, or lateral 
axis, is simply defined as being normal to the X and Z axes. Fortunately, 
due to the lateral symmetry, an airplane has no dynamic unbalance about the 
Y axis and thus only one product of inertia, I , exists. 

A second set of references, the earth reference system, must also be 
considered due to action of components of gravity along the three axes. 

To determine the gravity components, the airplane's attitudes in pitch and 
roll must be determined and the acceleration of gravity resolved along the 
three axes by utilizing the Euler equations. The pitch and roll angles are 
illustrated In figure 2, where pitch angle is show with wings level and 
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Z 
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1. Forces acting on the airplane 
along the axes: 
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2. Moments acting on the airplane 
about the axes: 
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ft- lb. 


3. ( Moments of Inertia of. the 
airplane about the axes: 


lx 


iy 
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* \ 

slug-ft^ 


4. Linear velocities of the 
airplane's CG in space: 


u 


V 


w 


ft/sec 


5. Angular displacements of the 
airplane about the axes: 


f" 


e 




rad. 


6. Angular velocities of the 
airplane about the axes:.. 


p 


q . 

. J 


r 


rad/sec 



Figure 1. Body axis system and parameters used in the 
digital simulation 



t 



3 



In these two row! it Ions the earth 



ro 1 { u< 1 1° sac #n i .. I ■(> *!y 'ev.^ 

reference angles a re. equivalent to Integrations of airplane angular vel- 
ocities about their respective body axes. But, if the craft is in an 
attitude with non-zero values of both pitch and roll, these angles must be 
obtained from the projection of the X and Y axes onto two vertical planes 
containing these axes, 

A third earth reference angle which has not been illustrated and which 
plays no part in the dynamics of motion is the yaw, or heading, angle which 
is the projection of the longitudinal body axis onto the horizontal plane. 

Characteristic Airplane Dynamics 

The six degrees of freedom of motion of an aircraft can now be defined 
as three angular freedoms in roll pitch and yaw about the X,Y* and Z axes 
respectively; and longitudinal, lateral, and normal translations along the 
X, Y, and Z axes respectively. Thus, the dynamics can be represented in 
terms of six simultaneous differential equations of motion. Before dis- 
cussing the makeup of the six equations of motion, it is felt that a better 
understanding of the following material may be obtained by a review of the 
various modes of aircraft motion. The motions can generally he divided in- 
to two classes, longitudinal modes and lateral-directional modes. The first 
class involves pitch rotation combined with translations along the X and Z 
axes, whereas the lateral-directional modes involve rolling* yawing^ and 
translation along the Y axis. This division of modes occurs because of the 
low level of aerodynamic and inertial cross coupling between longitudinal 
and lateral motions. For small amplitude oscillations (excursions in pitch, 
roll, and yaw less than 15-20 degrees) the two classes can be assumed dynami- 
cally independent . 
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Figure 2. Earth reference system used in the digital simulation 
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"1 -'*ii <2 . s • * c s m <>ro i1or« ■' fr 'ie ’’ta ,iiu’Lnal a*'d lateral 

6 irscCior.nl mo las of a rv resentutlve supersonic a Lrp lane <> The top half 
of the S plane is shown for each case. It is seen that the longitudinal 
motions occur in two oscillatory modes* one with a relatively long period, 
the phngoid mode, and one with a much shorter period, the short-period 
mode. The damping of the phugoic? is usually lot* and in some cases it is 
a neutral or even divergent oscillation. However, due to its long period, 
the motion is easily kept under pilot control with little or no ccxiscious 
attempt at correction, and no artificial stability augmentation is re- 
quired for this mode. The phugold occurs at constant angle of attack and 
manifests itself as a continuous trade-off between airspeed and alt tude. 

Figures 4. and 5. show that the period and damping ratio of the short 
period mode o£ a supersonic aircraft can vary over wide limits, and the in- 
clusion of the time lag of the human pilot in the control loop can some- 
times cause the short-period poles to move into the right hand side of the 
S plane. This necessitates the use of pitch stability augmentation not 
only for improvement of handling qualities but for reasons of safety. 

This mode occurs at constant airspeed and appears as a pitching oscilla- 
tion accompanied by small amplitude translations up and down the Z axis. 

It is extremely important in maneuvering flight and is of prime interest to 
designers of airborne fire control and stability augmentation systems. 

The lateral-directional pole* zero plot for a supersonic aircraft 
normally exhibits only one pair of complex poles and tvc real poles as 
shown in fig. 3. The frequency and damping ratio of the lateral-direction- 
al short period, or Dutch roll, mode vary as flight conditions change in a 
manner similar tc the variation of longitudinal short-period characteristics. 
This mode is a combined rolling-yawing oscillation with small amplitude 
lateral translations, and its name is derived from the swaying motions of 
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Figure 3 Pole-Zero Plot of the Characteristic Modes of 
Motion of a Supersonic Aircraft. 
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a Dutch skater. The lutch roll oscillations of a supersonic aircraft 
usually require artificial clamping in the same manner as required by the 
longitudinal short-period mode. Therefore, even though this study is con- 
cerned with an adaptive longitudinal stability augmentation system, the 
results and conclusions are applicable to the lateral-directional case 
as well. 

The roll damping, or roll subsidence, mode is shown as a negative 
real pole and it represents the viscous roll damping of the aircraft. 

And, the spiral mode is usually a steady divergence which manifests itself 
as a gradually steepening turn with accompanying sideslip. The divergence 
is fortunately slew and, as in the case of the longitudinal phugoid mode, 
is easily controlled by the pilot. The roll subsidence and spiral mode 
poles can occasionally form a complex pole pair in the case of some sub- 
sonic airplanes with large dihedral effect. This lateral-directional 
"phugoid” mode will exhibit small amplitude, long period rolling oscillations 
with considerable lateral translation and very little yawing oscillation. 

This concludes the survey of normal modes of airplane motion. How- 
ever, there are other "pathalogical" modes which can be experienced when 
the rates and amplitudes of motion become large. Some of these are the 
phenomena of inertial cross coupling in roll, autorotation (Ref. 125.)* 
and longitudinal pitch-up (Ref. 86); all of which can be eliminated or al- 
leviated by the use of well designed airplane stability augmentation systems. 
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Aerodynamic Stability Derivatives 
One final topic, the concept of aerodynamic stability derivatives, 
should be mentioned before proceeding to the equations of motion. The 
forces and moments which act on an aircraft are generally expressed in 
terms of non-dimensional coefficients multiplied hv appropriate aircraft 
and flight condition parameters. As an example, the equation for airplane 
lift is shown below, 

P 2 

Lift * C L j VS 

Where : 

Lift is in pounds 

C is the airplane lift coefficient 

Li 

p s atomspheric density in slugs/ft~ 

V 55 airplane forward velocity in ft/sec e 

2 

S = airplane wing area in ft . 

P 2 

The expression — V is a measure of the kinetic energy of the air- 
stream and is termed the dynamic pressure because it also represents the 
force which would theoretically be exerted against a unit square surface 
held normal to the stream. Thus, the product of dynamic pressure and wing 
area represents a hypothetical aerodynamic force from which the actual air- 
plane lift can be obtained through multiplication by the lift coefficient. 
Moments acting about the aircraft can be obtained in the same manner with 
the inclusion of a moment arm, for dimension considerations, as is done in 
the following pitching moment equation. 

P 2 — 

Pitching Moment » C t V Sc 
° m2 

Where the moment arm c is the wing mean aerodynamic chord. Corres- 
ponding expressions for drag force, side force, rolling moment, and yaw- 
ing moment coefficients are included in the table of symbols. It is seen 
that there exists an aerodynamic coefficient for each of the six degrees 
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of freedom or motion, These coefficients and their derivatives are 
usually determined empirically in wind tunnel and flight studies be- 
cause of the inaccuracy of analytical determinations. 

Each coefficient can be expressed in a Taylor f s series such as: 



Cr m 

•••+(&) 



<Z + \dCC 



M Ct.\ ao 1 

\dCC L l 2.1 



■-H 



(c )C L \ . 

\7gTioc+ ladTaj 



oc ‘ 

\dcc*J 2f 



+ ■ 



2 ! 



+ 



The effects of the second order and higher terras are insignificant in aero- 
dynamics and are dropped. So, for practical purposes the coefficient can 



OC 



4- 



IoCl 






be exnressed: Ic)Cl\ , OCl\ 

Cl=1 3 <xI oc + ( a< xj 

The partial derivative terms such as are called aerodynamic stability 

derivatives and are usually written for brevity in the formC^^ 

• . C L = C Lcc OO d- 6 lo> CL + C iq 

Figure 6, a curve of lift coefficient vs, angle of attack, can be used to 

illustrate the physical meaning of stability derivatives. Here,( ~ ^r- ) 

is the slope, or derivative, of the airplane lift curve. It is seen that 

the curve is linear up to <X ^ 12°. This indicates that AC ~ C CC is 

L hoc 

physically valid for small variations in OC » The fact that almost all 
aerodynamic stability derivatives remain nearly linear over a similar range 
of their corresponding variables enables the small amplitude motions of the 
aircraft to be expressed as a set of simultaneous linear differential equa- 
tions. 

How many stability derivatives must be considered in the equations of 

k k 

motion 9 Given the six force and moment coefficients , C xr> , Ch , 



*The lift coefficient and the drag coefficient act perpendi- 
cular and parallel, respectively, to the relative wind, i>ut it can be 
assumed with little error that they are also normal and parallel to the 

X body axis. Then C_ = C„ and C = €_« 

Z u a D 
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Angle of Attack, deg. 

Figure 6. Plot of Airplane Lift Coefficient vs. Angle of Attack. 

I 

and C and the fact that each coefficient can be "inf luenced by: 



1. angular displacement about each axis, 

2. angular rate about each axis, 

3. velocity along each axis, and 

4. acceleration along each axis 



there are 12 possible derivatives for each coefficient or 72 in all. Also, 



control surface deflections can add 36 more terms. ^'Fortunately, airplane 
symmetry, which causes a lack of cross coupling between longitudinal and 
lateral - directional modes, greatly reduces the number. The magnitudes of 
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some of the remaining derivatives are so small as to make them insigni- 
ficant, As a result there are approximately 25 derivatives considered 
important enough to he included in a simulation program. These are sum- 
marized in Appendix TIT, A list of stability derivatives for a representa- 
tive supersonic aircraft is tabulated in Appendix I. The aircraft is the 
same one for which the longitudinal short period mode characteristics are 
illustrated in figures 4, and 5, 

The Airplane Equations of Motion and the Longitudinal Short Period Transfer 
Function, 

The six equations of motion corresponding to the six degrees of free- 
dom of an aircraft are derived in appendix IV and the final results are 
presented in table I, It should he noted that the equations represent the 
balance of inertial and aerodynamic forces and moments along and about the 
airplane body axes, and that the terms are composed of .airplane physical 
constants, aerodynamic stability derivatives, and aerodynamic variables. 
Assumptions used in the derivation of the equations are: 

1. The change in aerodynamic forces and moments due to deviations 
from an equilibrium condition are assumed directly proportion- 
al to the deviation. 

2. The deviations of the airplane from equilibrium are assumed 
small, 

3. The airplane is a rigid body. 

4. Constant thrust is assumed during an airplane distrubance, 
and secondary power effects are neglected. 

The equations are used directly in the digital computer simulation of the 
airplane dynamics, which is outlined in appendix VI S and the airplane 
longitudinal short-period transfer function is derived from the pitching 
and normal translation equations of motion in appendix V. This transfer 



function is of prime importance in the study of the airplane pitch damper 
because the short -period mode is the dominant longitudinal mode. The 
short-period transfer function is written below both in its general form and 
also in terms of the airplane physical characteristics, airplane stability 
derivatives, ar.d aerodynamic parameters. The second form of the transfer 
function is shown to illustrate the dependence of the dynamics of the short 
period mode on changes in stability derivatives and aerodynamic parameters, 
and it is possible from inspection of the function to determine which de- 
rivatives and parameters exert influences on a given characteristic, such 
as damping ratio, of the short-period mode. 



Longitudinal Short Period Transfer 



Function 



A. General Form 
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B. Detailed Form 
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sin^> cos 0 & pitch angle 

5 o roll angle 

cos® cos d , _ _ . 

' (Angles between the body axis 

and earth gravity axis) 
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Ill - THE !ttB OF A FITCH l ' : 0 L0L40 IT UT) I hVJ. ST. mm V S •LifTATIOH 

The purpose of airplane stability augmentation is to reluce the 
range of variation of a dynamic characteristic, or characteristics* and 
to hold the variables of interest ns close as possible to predetermined 
optimum values* In the case of an aircraft these optimum values are de- 
fined in terms of "airplane handling qualities'* which describe an air- 
plane* s response to manual control by a pilot. Even though the determina- 
tion of desirable handling qualities is subjective, and is dependent on 
pilot opinion, a number of studies have been conducted to specify these 
qualities in quantitative engineering terms* One of the better known 
studies, conducted at the Cornell Aeronautical Laboratory, (Ref* 34), con- 
cludes that the desirable natural frequency and damping ratio of the short- 
period mode of a fighter airplane lie in the range CO n s .4 to *5 cycles per 
second and damping ratio ~ .65 to ,70. The same study also concludes that 
even though the longitudinal natural frequency is of importance in maneuver- 
ing flight, it can vary over wide limits and the handling qualities will 
still be considered tolerable, as long as the short period damping ratio is 
held in the .65 to .70 range* Thus it is considered more desirable to equip 
a high performance airplane with artificial damping than to provide static 
stability augmentation. It should be noted that for the craft considered 
in this study, as well as for almost all other airplanes of similar type 
and mission, the longitudinal damping ratio of the un-augmented airplane 
is always below the desired .7, as is shown in figure 4. This means 
that the pitch damping system will work to supplement the basic aerodynamic 
clamping in all cases. 
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Mechanization of the I itch hamper 



In the discussion of the airplane longitudinal short “period transfer 
function it was shown that the short-period motion can be represented as a 
second order system with characteristic damping ratio and natural frequency. 
It was also shown that the natural frequency and damping ratio can be ex- 
pressed as functions of the aerodynamic stability derivatives. The ex- 
pression for short period damping ratio being: 



From the above relationship it can be seen that the most effective means 
of varying the damping ratio is to artificially vary the pitch damping de- 



resistance to pitching velocity. 

The pitch damper must then provide an aerodynamic pitching moment which 
is proportional to airplane pitch rate, and which has a direction in opposi- 
tion to the angular rate. Thus, the moment produced is in phase with the 
viscous aerodynamic damping in pitch and the two moments combine to form 
the total pitch damping moment. The artificial damping moment is usually 
generated by deflecting the airplane 0 s horizontal stabilizer in direct 
proportion to measured pitch rate. 

Figure 7. is a block diagram of the pitch dampen u r u:A In this study* 

The airplane transfer function shown in the diagram represents the longi- 
tudinal short -period mode only because, as will be demonstrated later, the 
phugoid mode and the lateral-directional dynamics are of negligible impor- 
tance in longitudinal stability augmentation considerations. 

The stabilizer actuator in the diagram is an irreversible hydraulic position 
servo, and the parameters in the transfer function represent the state~o£- 
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Figure 7« Block Diagram of the Pitch D&mpar 
the-art. as used in modern aircraft. It is generally assumed that the 
dynamics of the rate gyre have negligible effect on system operation be- 
cause the ratio of the characteristic frequency of state-of-the-art rate 
gyros to the short-period mode frequency is of the order of 20 to 30 
(Ref. 43), and thus the gyro transfer function is assumed to he unify. 

Table II lists the II airplane transfer functions use I in the study alon^ 
with corresponding flight conditions. And, figures S thru 18 are root 
locus plots showing the behavior of the system with varying dt-mp^r gain at 
each condition. The square symbols on the plots indicate the gair values 
at which the airplane short-period damping ratio has been increased to 0.7. 
It is seen that the use of the carper has %* 2 ty little influence on the 
short period frequency when operate \ over the range of gains required to 
bring £ up to 0.7 at each flight condition* And, it is also seen that 



the actuator 



pole moves i r, such a direction, with increasing damper gain, 
that the actuator damping ratio regains close to 0, 7 over the danper opera- 
ting range* The actuator rode will become unstable at extremely high gains* 
but operation at such gains is of no interest in this study* 

It is known that the closing of a feedback loop around a plant having 
varying characteristics will cause the plant output to vary over a smaller 
range, as plant characteristics vary, than the range of variation of out- 
put with no feedback* From this it would be expected that the pitch damper 
should not only Increase the average value of damping ratio over all flight 
conditions, but also that the range of variation of damping ratio, as condi- 
tions vary, should decrease* This is borne out in Figure 19. which is a 
plot of short-period damping ratio vs, speed and altitude with constant 
damper gain s 0.3. The chosen value of damper gain is approximately the 
median of the gains required for 0*7 damping ratio over all 11 flight 
conditions. In Che case of the basic airplane, £ varies from 0.445 at M 
= 0,9 <? 10,000 ft. to £ = 0.064 at M = 1*2? 45,000 ft. whereas the varia- 
tion of / vrith damper gain » 0.3 is from 1,1 to 0.436 for the same two 
flight conditions* Thus the range of variation about the mean has been 
reduced from approximately 190% to approximately 95%. 

Even though the use of the constant gain damper has caused a consider- 
able improvement in the dynamics of the airplane, a further restriction in 
the range of variation of damping ratio is usually desired for fighter-inter™ 
ceptor type aircraft* One method of solving this problem, which is in 
general use today, is the programming of stability augmentation system 
gains, as functions of measured flight parameters, in such a manner that 
the programmed gains will cause near optimum system response over all 
flight conditions. The techniques of gain scheduling methods range from 
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variation oi Jasper gain nr function of e.casu red airspeed to systems 
which program as functions of airspeed* altitude, fuel load, center of 
gravity position, wing flap position, landing gear position, and vari- 
ous combinations thereof. One fact becomes obvious from the preceeding 
discussion; the dynamic characteristics of the aircraft must be accur- 
ately known in advance at all flight conditions and for all airplane 
configurations if the programming system is to operate successfully * 

Such a detailed determination of airplane characteristics requires a 
large amount of wind tunnel and flight testing which is both time consuming 
and costly* And, even after the desired information has been determined 
it is usually necessary to conduct lengthy flight tests with the program- 
ming and stability augmentation systems installed in the airplane for final 
"tailoring*' of the aircra ft- system combination* In addition, accurate in- 
flight air data and airplane configuration Information must be continusouly 
fed into the system during flight to provide the basis for gain scheduling* 
The instrumentation and transducers required are generally costly, and are 
usually complex enough that they create reliability and fail-safety problems, 
A possible solution to these control problems is to design a system 
which is intelligent enough to determine when the airplane °s dynamic per- 
formance has changed, and in what manner it has changed, and to use this 
information to adjust the stability augmentation system gains to obtain 
optimum performance over all flight conditions* This is the concept of 
the self-adaptive control system. 
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IV 



SELF-ADAPTIVE CONTROL SYSTEMS 



The concept of adaptive control has been applied to a wide range of 
control system types with the result that a number of definitions of adap- 
tive control exist. However, the most often used definition* and also the 
one which appears to best fit the aircraft control problem is as follows: 

An adaptive system is a control system configuration in which 
the measurement of process dynamics or signal characteristics 
is utilized to automatically adjust the controller in an attempt 
to achieve optimum operation at all times. 

From this definition it can be seen that the adaptive process must per- 
form two functions, the first of which is the ,f identification” function, 
or the measurement of process dynamics cr signal characteristics to deter- 
mine present system state and compare it with the desired state. And, the 
second is the "command” function which is automatic changing of the system 
so that it will approach a desired state. In most adaptive control studies 
it is assumed that the variation of system characteristics takes place slov; 
ly with respect to the frequency of the system mode of interest, and thus 
the identification process can be allowed a relatively long measurement in- 
terval during which to determine system state. However, in the ail craft 
problem, it is desirable to perform the adaptation process relatively rapid 
ly in comparison to the quantity being controlled, in this case the short 
period mode dynamics. This is because virtually instantaneous changes can 
take place in the airframe transfer function due to release of payload or 
certain Mach number effects on stability derivatives. For this reason* 
one of the objectives of this study is the investigation of an adaptive 
technique capable of rapid identification and command. Another considera- 
tion in mechanization of the adaptive process is the generally accepted 



requirement feu* 3 > iu o :Iei :oi the system it be con.si.deve 1 truly adap- 
tive, the identif itaiion process can net involve measurement of any 
variable which Is external to the system. To illustrate the concepts 
and requirements discussed above, a block diagram of Gibson* s ideal 
generalised adaptive system (Ref. 76) is shown in Figure 20 . 
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Figure 20, Gibson 1 s Ideal Generalized vlaptive 

System 

Classification of Adaptive Systems 
In o'cdei' to compare the approaches which have been taken to the 
design of adaptive control systems, Aseltlne, in reference 4* has seperat 
ed them into 5 classes which, in some cases, have a certain amount of 
overlapping, but even so, they form a convenient framework for comparison 
Aseltine's classes are: 
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1) assive Vi i tat: ion: Fastens winch adapt without automatic 

changes ir» system characteristics, but rather through de~ 
sign for uniform operation over wide variations in environ- 
ment . 

2) Input Signal Adaptation: Systems which automatically ad- 

just their parameters in accordance with input signal 
characteristics* 

3) Sxtrenam Adaptation; Systems which adjust to obtain a maxi- 
mum or minimum of some system variable, 

4) System Variable Adaptation: Systems in which the command 

process is actuated directly by measurement; of a system 
variable or variables* 

5) System Characteristic Adaptation: Systems in which the 

identification process measures system variables in order 
to determine system characteristics. The system character- 
istics thus determined are compared with desired character- 
istics and the differences used to generate command signals* 

It has been noted from a survey of adaptive systems which have been invest 
gated for possible use in aircraft control that, to date, all such systems 
fit into classes 1, 4, or 5. And, it has also been noted that an alter* 
nate classification system might be applied to aircraft adaptive control 
systems. The suggested classes are: 

1) High Gain Adaptation: System s in which the characteristics 

of the aircraft are submerged by the use of high feedback 
and/or forward loop gains. Some form of idealized model is 
normally used to provide either system input signals or out- 
put comparison signals, and the techniques usually involve 
either nonlinear processes, with attendant limit cycle, or 



operation of a minor system mode near instability. 

2) Impulse Response Adaptation; Systems in which the identi- 
fication process computes a time history o£ the plant im- 
pulse response. The impulse response is then utilized by 
the command process to generate plant alteration signals, 

3) Parameter Adaptation: Systems in which plant parameters 

are determined by direct computation. The error between 
the desired value of the parameter and its actual value is 
used as the command signal for plant alteration. 

Survey of Adaptive Techniques which have Been Applied to the Airplane 

Stability Augmentation Problem 

1) High Gain Adaptation Systems 

The earliest studies of aircraft adaptive control involved systems 
which operated at high gains in order to force the variation of airframe 
parameters into a small range. One of the first of these studies was 
conducted by Rath of WADC (Ref. 134) on a method which was not truly ad- 
aptive, in that parameters were not automatically varied* however* it 
may be classed as passive adaptation. Rath* s, system used a parallel com- 
plementary optimum response model as shown in Figure 21. The goal of 
this study was to design a system which would produce optimum airframe 
dynamic characteristics, as airframe parameters varied, with constant 
system gain. A fourth order transfer function was selected as the opti- 
mum model for the pitch dynamics and the same input signal was applied to 
the airframe and to the model. The difference between model and plant 
output was obtained and various functions of the difference signal were 
utilized to improve the airframe response, Both the forward loop gain 
and the shaping loop gain were made high in order to M force” the system 




Figure 21. Rath*s System for passive Adaptation Using a Complementary 

Optimum Response Model 

overall characteristics to resemble the model. Feedback functions studied 
included pitch angle* or position* error feedback, integrated position 
error feedback* rate error feedback* and rate-error-plus-position-error 
feedback. It should be noted that error* as referred to here* means the 
difference between plant output and model output rather than plant input- 
output error. The first two feedback functions produced overall instabil- 
ity at high gains* whereas rate error feedback resulted in a stable system 
with infinite steady state position error. Rate-error-plus-position-error 
feedback allowed very high gains to be used before the system became un- 
stable* and while the system was stable* the airframe output rr&tcbed the 
model output acceptably at most flight conditions while using a constant 
value of feedback gain. It was also found that* whatever feedback shaping 
was used, the basic airframe and hydraulic actuator must fom a stable 
system with the shaped feedback loop alone. That is* the feedforward loop 




could rot coft cc m-st * >L {■ > of rhe a.-.’ shaped fV^Jhack 

combination. The optimum **ra *?fer function feedforwrd in combination 
with the shape! fee-ib^c 7, h**d the. effect of changing the input to the 
airfrai a to a v&Iua which caused 'the airframe to have the desired out- 
put, assuming the airframe htl th* desired control rates and high gain 
control power « 

An approach similar ro RairhAs passive adaptation Method was taken by 
Jandols of Corwair who utilised the inverse of the optimum response model 
in the feedback loop in combination with a high forward loop gain system 
as shown in Figure 22* (Ref* 43). 



AIRFRAME SHORT PERIOD 
TRANSFER FUNCTION 




Figure 22. Dandois 8 System Using the Inverse of the Optimum Response Model 

in the Feedback Loop 

He demonstrated that, with such a feedback arrangement, the overall system 
characteristics approached the characteristics of the optimum response 
model as forward loop gain was increased. In actual operation, the 
system held air r rana parameters wichir accet table lir.its c»&r - wide 



range oC flight con ilt ions and the i\ghea. the gain the better the 
performance, However, gain 1 units, dependent upon "light conditions, 
were reached at which the systen became unstable. The instability was 
shown to be due to sensor non linearities rather than a basic, theoretical, 
system instability. 

One of the first "high gain" adaptive systems to actually use an 
identification and command process was the system-variable-adaptation 
system of Fluggc-Lotz and Taylor (Ref, 64), which is Illustrated in Fig* 
ure 23, As is shown in the block diagram, one of a number of alternate 




Figure 23, Nonlinear Adaptive System of Flugge-Lotz and Taylor 
feedback paths was used according to a switching logic process based on 
simultaneous measurements of system output, error signal, and derivatives 
of these. The Identification process thus consisted of the measurements 
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*n> sv ■ tch* "a : «. ••*» u • rb • *.pj. rvf u/tb ^e«r luaU path, and 

the coniP-.ati^ cons ■ sle of in-, itching to the chosen path. Forward gain was 
kept high in order to “force” a small instantaneous error between input 
and output, and the system was nonlinear during any length of time com-* 
parable to the natural period of the plant because of the switching be- 
tween feedback loops during one cycle of operation* The arrangement was 
shown to be insensitive to variations in plant parameters over wide limits, 
however a limit cycle persisted during operation due to the feedback non- 
linearity. 

The Minneapclis-Honeyvel! relay actuated adaptive system (Ref* 165), 
which is presently undergoing flight tests, is similar to the Flu&ge-hotz 
mechanization but uses only one feedback path with switching logic in the 
forward loop. The switching logic derives Its input from the error signal 
between optimum response model output and the process feedback signal as 
shown in Figure 24* The switching logic actuates a relay, or on-off con 
troller, to operate the process at a high gain with the objective of 
maintaining a small error signal between model output and plant output* 

As in the Flugge-Lots system, the airframe output approximately dupli- 
cates the model output when the forward gain is high, and the character- 
istic limit cycle is also present. A second adaptive feature of this 
system, aside from the switching logic, is the us^ of limit cycle ampli- 
tude measurement to vary forward loop gain. That is, forward gain is ad- 
justed in such a manner that the limit cycle amplitude is held constant 
over all flight conditions. The prototype of this system was successfully 
tested in an F-94C fighter aircraft and an improved version is presently 
being flight tested in the X-15 high altitude research airplane. 

Another application of the utilization of the inverse optimum re- 
sponse transfer function ixi the feedback loop, in conjunction with 
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Figure 24. MJnneapolis-Moneywe! 1 Relay Actuated Adaptive System 
high forward loop gain, is the General Electric system (Ref. 32), which 
was recently tested in an F-102 interceptor aircraft. The mechanization 
is similar to that of Dandois shovn in Figure 22, except that ar> active 
system-variable-adaptat ion feature has been added. This feature Involves 
variation of the forward loop gain, as flight conditions change* to 
maintain the system just below the boundary of instability. As was pro- 
viously mentioned, when considering systems of this type* the higher the 
gain, the closer will be the airframe characteristics to the optimum re- 
sponse characteristics and thus if the system Is operated near the in- 
stability boundary, the maximum performance will he extracted. In opera- 
tion, the adaptation ^ stem monitors the frequency the control surface 
actuator because the actuator mode is the one which becomes unstable as 
forward gain is increased. The g*in Is then adj isted to iraintn Ic the 
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of this rode ecus 



It should 



frequency of this mode constant over a 1 1 flight conditions, 
be noted that: the actuator mode is a lightly damped oscillation and 
that normal control Inputs are utilized to provide the required pertur- 
bations. The frequency measurement is accomplished with a narrow band 
pass filter and zero crossing counter arrangement. 



2} Impulse Response Adaptation Systems 

The function of the identification process In the impulse response 
adaptation class of systems is first, the determination of a time history 
of the plant impulse response through measurement of system input and out- 
put signals over a period of time, and second, the determination of plant 
characteristics through examination of the impulse response. An early 
method of impulse response self optimization which was investigated by 
Aseltine (Ref. 4} serves to illustrate the process, even though this 
particular mechanization would not be suitable for aeronautical use. 

The impulse response was obtained by periodic input pulses which were 
widely spaced with respect to the period of the plant mode of interest. 
After each perturbation, the system damping ratio was computed by separ- 
ate integration of the positive and negative portions of the inpulse re- 
sponse and obtaining the ratio of positive to negative area. Damping 
ratio was computed from ibis ratio and the system damping coefficient was 
varied in a closed loop manner to cause actual damping ratio to approach 
a desired value. This met he 1 was considered undesirable for aircraft use 
primarily because normal operating signals must be excluded from the in- 
put during the comparatively lengthy identif icatior period. 

A method of impulse response determination in which normal Input 
signals could be included In the identification process was investigated 



by / nee r son c «: h.ivj? -in’ s ^ ’ » fho *ech*vii.?ai; ion «ane use of the 

convolution theorem wheuzby the system input convolved In t ini'; with the 
impulse response lesuHs in a time history of the system output* In opera- 
tion, random noise 'with auto*: orrel sti cm Z*mrtto\\(pii('Y) was superimposed on 
the norma! input signal is shown in Figure 25. The plant output was filt- 
ered to separate the noise component front* the tot a J resj o>nse» and the in- 
put and output noise ccrx.t. uents were cross-correlated to obtain tf)jo (T) - 

CO 

. The noise input had a bandwidth 3 to 10 times the pj^rit 

co 

bandwidth so that ^ f)ii(T ') was af fectively impulse and (pio(/) an - 

proximate ly the plant impulse response. Twelve channels of corss cor re la- 
lion were used and e^ch channel provided one point on a time history of 
the impulse response. Damping ratio was determined an:I utilized as its - 
cussed in the proceeding paragraph. 
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Figure 25. 



Block Diagram cf :ercontronic*s Impulse Ttesporse c e]f {Vtimi- - 

i vg System 
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